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This  paper  applies  the  spheroidal  T-matrix  approach  of  Dr.  Roger  H  Hackman  described 
in  his  paper  “The  Transition  Matrix  for  Acoustic  and  Elastic  Wave  Scattering  in  Prolate 
Spheroidal  Coordinates”,  [Journal  of  the  Acoustic  Society  of  America,  75(1),  PP.  35-45]  to 
describe  the  scattering  from  large  aspect  ratio  targets.  Due  to  the  limits  of  double  precision 
arithmetic  in  the  computation  of  the  spheroidal  wave  functions,  these  methods  are  limited  to 
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I.  INTRODUCTION 


Prolate  spheroidal  coordinates  are  one  of  eleven  coordinate  systems  in  which  the  scalar 
Helmholtz  Equation  in  three  dimensions  is  separable.  Separability  is  the  reason  the  scattering 
from  a  rigid  (sound  soft)  prolate  spheroid  may  be  expressed  as  a  partial  wave  series  of  spheroidal 
wave  functions  analogous  to  the  partial  wave  series  of  a  rigid  sphere.  However,  the  vector 
Helmholtz  Equation  is  not  separable  in  the  prolate  spheroidal  coordinate  system.  This  is  the 
main  reason  for  the  lack  of  exact  solutions  for  acoustic  scattering  from  an  elastic  prolate 
spheroid.  The  lack  of  separability  of  the  vector  wave  equation  in  spheroidal  coordinates  is  also 
the  reason  there  does  not  exist  an  expansion  of  the  Green’s  Dyadic  for  the  vector  wave  equation 
in  terms  of  spheroidal  functions. 

In  the  early  1980’s,  Roger  Hackman  1-3  at  Coastal  Systems  Station  used  Betti’s  Identity  to 
derive  a  T-matrix  description  of  the  scattering  from  large  aspect  ratio  elastic  targets  based  on  the 
spheroidal  wave  function,  rather  than  the  popular  spherical  T-matrix  description.  This  approach 
is  numerically  more  stable  than  the  spherical  T-matrix  approach  for  large  aspect  ratio  targets. 

The  author  had  the  pleasure  of  working  with  Dr.  Hackman  during  his  tenure  at  COASTSYSTA, 

A  Q 

and  worked  on  the  scattering  from  large  aspect  ratio  targets  in  a  waveguide  '  . 

This  article  focuses  on  the  use  of  the  spheroidal  T-matrix  description  of  scattering  from 
elastic  targets  to  describe  the  low  frequency  scattering  from  these  targets.  The  purpose  for  this 
work  is  to  build  an  object  oriented  computer  program  for  computing  the  low  frequency  scattering 
from  large  aspect  ratio  targets.  This  computer  code  is  to  be  integrated  into  the  sonar  simulation 
PC  Shallow  Water  Acoustic  Tool-set  (PC  SWAT)  to  describe  a  low  frequency  imaging  sonar  in 
the  presence  of  elastic  returns  from  the  target. 

In  addition  to  the  Introduction  and  Reference  sections  of  this  report,  the  outline  of  this 
paper  includes  the  following  sections.  Section  II  describes  the  spheroidal  coordinate  system  and 
its  relationship  to  linear  acoustics.  Section  III  outlines  the  computation  of  the  spheroidal  basis 
functions.  Section  IV  describes  the  computation  of  the  spheroidal  T-matrix.  Section  V  contains 
a  collection  of  sample  calculations  and  a  comparison  of  the  author’s  results  with  those  of 
Hackman.  Most  of  the  details  of  the  computations  are  contained  in  the  appendices.  Appendix  A 
describes  the  computation  of  the  eigenvalues  and  expansion  coefficients  for  the  spheroidal  wave 
equation.  Appendix  B  contains  the  orthogonality  and  completeness  relationships  of  the 
expansion  coefficients.  Appendix  C  contains  the  expansion  of  the  scalar  Green’s  Function  and 
plane  wave  in  terms  of  spheroidal  wave  functions.  Appendix  D  describes  the  transformation 
between  spherical  and  spheroidal  wave  functions.  Appendix  E  describes  the  computation  of  the 
radial  spheroidal  functions.  Appendix  F  contains  a  tabulation  of  the  connection  components  and 
its  first  order  derivatives  used  in  the  computation  of  the  second  and  third  order  covariant 
derivatives  of  the 
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scalar  basis  functions.  Appendix  G  contains  a  tabulation  of  the  components  of  the  vector  basis 
functions  and  their  covariant  derivatives.  Appendix  H  contains  a  tabulation  of  the  parity  of  the 
basis  functions.  Appendix  I  describes  the  computation  of  the  surface  integrals  and  T-matrix  in 
the  case  of  a  prolate  spheroid.  Appendix  J  describes  the  computation  of  the  surface  integrals  and 
T-matrix,  in  the  case  of  a  finite  cylinder  with  hemi-spherical  end  caps. 


2 


CSS/TR-02/18 


II.  SPHEROIDAL  COORDINATES 


Prolate  spheroidal  coordinates  are  a  generalization  of  the  familiar  spherical  coordinates. 
The  transformation  from  spheroidal  to  Cartesian  coordinates  is  given  by  the  following  equations. 

*  =  /V(f’— W-'rtcosW 

•V^AAf-Dd-^sinW  (1) 

2  =  /^ 

Here,  £  >  1  is  the  radial  coordinate,  -1  <  77  <  +1  is  the  angular  coordinate,  0<(p<lK  is  the 
azimuthal  angle  coordinate,  and /is  the  semi-focal  distance. 

Adopt  a  right-handed  coordinate  system  whose  coordinates  y*  are  given  by  the 
spheroidal  coordinates  (£ ,(p,7 ]).  The  metric  tensor  in  this  coordinate  system  is  a  diagonal  metric 
whose  diagonal  members  are  given  by  the  following  expressions  in  terms  of  the  coordinate  basis 
0|,dp,d9}  of  the  tangent  space: 


gs=/!(#! -tfVtf -l) 

(2) 

Define  the  ortho-normal  triad  {e,,e2,e3}  =  {%,((),  fj}  by  the  following  relations: 

s,=<3, 

ef^Sflh,  (3) 

hfi  ~  4J&w 


The  covariant  derivative  of  a  tensor  of  type  (1,1)  in  terms  of  the  coordinate  basis  is  of  the 
following  form: 
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Here,  V* ^  are  the  components  of  the  torsion-free  affine  connection  defined  by  the  following 

expression  in  terms  of  the  metric  tensor  (An  affine  transform  changes  the  relationship  between 
the  coordinate  system  a  program  uses  to  draw  and  the  coordinate  system  used  to  display). 

rV  =^ga<T(df<g,<T  +dvgM< y  —  ^agfiv  )  (5) 

This  connection  has  the  property  that  the  covariant  derivative  of  the  metric  tensor  vanishes. 

The  stress  tensor  of  an  elastic  solid  is  defined  by  the  following  expression  in  terms  of  the 
metric  tensor  and  the  covariant  derivative  of  the  displacement  vector  u. 

Tap  =  X  {yaua)gap  +M(VaUp  +Vp“a)  (6) 

Here,  X  and  fi  are  the  Lame  constants  of  an  isotropic  elastic  solid,  which  are  related  to  the 
longitudinal  and  shear  wave  speeds  by  the  following  relations,  where  p  is  the  density  of  the 
solid. 

X  =  p(v2L-2v2T) 

*)  V  ' 

M=pv; 

The  equation  of  motion  for  an  isotropic  elastic  solid  is  given  by  the  following  expression, 
where  T(u)afi  is  the  stress  tensor  for  the  displacement  vector  ua . 

(8) 

A  generalization  of  Green’s  Theorem  gives  rise  to  Betti’s  Identity,  which  states  that  given 
a  pair  ( u,v )  of  regular  solutions  of  Equation  8,  the  following  surface  integral  over  the  boundary  of 
the  enclosed  volume  ( V)  vanishes. 


(Qds{u  •t(v)-t(u)»v}  =  0  (9) 

(IV 


Here,  the  terms 

t(u)^T{u)Mana  (10) 

are  the  components  of  the  traction  of  the  displacement  u. 
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III.  SPHEROIDAL  FUNCTIONS 


The  scalar  Helmholtz  Equation  in  spheroidal  coordinates  is  of  the  following  form: 


(v„ v” + e  )V = g"  Via, + efg)V = 


1 


f) 

,  (.Z2-v2) 


(f-l)(l -i72) 


{a?.(f -i)a^+a,(i-772)5,^ 

a^+c2(^2-^2)^}=o 


(ii) 


g  =  Det(g/lv) 

Here,  c  is  the  dimensionless  product  of  the  wavenumber  k  and  the  semi-focal  distance  f  This 
wave  equation  is  separable  into  the  product  of  a  solution  of  the  radial  equation: 

3{  (f  - 1  $(je„,  -  a.,  (c)-c2?+  =»  02) 

the  angular  equation: 

=0  (13) 

and  the  azimuthal  equation.  The  spheroidal  coordinates  have  the  rather  interesting  property  that 
the  radial  and  angular  equations  are  isomorphic  under  the  interchange  of  the  coordinates  £  <->  ij . 
This  isomorphism  implies  a  duality  between  the  solutions  of  the  radial  and  angular  functions, 
that  is,  the  radial  functions  may  be  regarded  as  the  analytic  continuation  of  the  angular  functions 
to  the  domain  +1  <  rj  <  +°° .  In  particular,  the  radial  functions: 

Kml 

ne„(c^)=-Ls'«(c,^ 

Kml 


(14) 
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are  proportional  to  the  corresponding  angular  functions.  This  transformation  is  an  isometry  of 
the  metric  tensor,  that  is,  the  metric  tensor  is  invariant  under  this  transformation.  The  classic 
method  of  solving  the  angular  wave  equation  is  to  expand  the  angular  functions  of  the  first  kind 
in  terms  of  associated  Legendre  polynomials: 


•?>.*;)=  I  <(<0 C,0?)  <15> 

*=0,1 

Here,  the  sum  is  over  either  even  or  odd  indices,  depending  upon  whether  the  difference  ( l-m )  is 
even  or  odd.  The  series  expansion  given  in  Equation  1 5  is  a  solution  of  the  angular  wave 
equation  provided  the  expansion  coefficients  satisfies  the  following  recurrence  relations: 

«,2+(Am-4/Hm/+«W2=o 


m  _  (2m  +  *  +  2)(2/w  +  k  +  l)c2 
k  (2m  +  2k  +  3){2m  +  2k  +  5) 

J3?  =(m  +  k)(m  +  k  + 1)  O6) 

2  2 (m  +  k)(m  +  *  + 1)  -  2 m2  - 1 
+C  (2m  +  2*  -  l)(2/n  +  2*  +  3) 


*(*-  Qc2 

Jk  (2m  +  2k-  3)(2m +2k-l) 

The  eigenvalue  A m,(c)  is  chosen  such  that  the  limit 


(<M 

d?(c ) 


)  =  0 


(17) 


vanishes.  The  conventional  method  of  solving  for  the  eigenvalue  is  to  solve  a  transcendental 
equation  defined  in  terms  of  a  continued  fraction  of  the  above  eigenvalue  and  coefficients  as 
described  by  Flammer  in  Reference  9.  A  numerically  more  robust  method  of  solving  for  the 
eigenvalues  is  to  transform  the  above  recurrence  relation  into  the  problem  of  solving  for  the 
eigenvalues  of  an  infinite  dimensional  tri-diagonal  matrix.  This  method  is  due  to  Hodge  l0,  and 
is  described  in  References  10  through  12.  This  is  the  method  adopted  by  the  author  since  it  is 
readily  generalized  to  the  case  of  complex  wavenumber.  The  coefficients  df  (c)  are  related  to 

the  eigenvectors  of  this  infinite  dimensional  matrix.  Reference  1 3  contains  a  simple  iterative 
method  for  solving  for  the  eigenvalues  and  eigenvectors  of  arbitrarily  large  tri-diagonal  matrices. 


In  the  case  of  the  radial  wave  equation,  the  fact  that  the  radial  functions  are  proportional 
to  the  analytic  continuation  of  the  angular  functions  can  be  used  to  calculate  the  radial  functions 
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for  radial  coordinates  near  unity  (+1  <£  <1.5)  as  described  in  Reference  1 1 .  The  computation  of 
the  irregular  radial  functions  requires  the  computation  of  the  angular  functions  of  the  second, 
which  is  derived  by  taking  the  following  limit  of  the  series  expansion  in  terms  of  associated 
Legendre  polynomials  of  the  second  kind. 

Cfo<7)  =  linv*  £  <,(<0fiL,„(!7)  (18) 


The  associated  Legendre  polynomials  of  the  second  kind  exhibit  a  simple  pole  at  integral 
orders  less  than  -m.  At  the  same  time,  the  coefficients  df(c)  exhibit  a  simple  zero,  and  the 
product  remains  finite. 

11  k>+2m  (19) 


This  produces  a  representation  of  the  angular  function  of  the  second  kind  that  is  a  series 
of  the  associated  Legendre  functions  of  the  first  and  second  kind  of  the  following  form: 

S%(c,r,)=  X  ftCW*  Z  (20) 

k=~2m+S  k=2m+2-S 


The  coefficients  d^p  (c)  are  defined  by  the  limit: 


<(c)  =  lim,_0 


(21) 


A  recurrence  relation  for  these  coefficients  is  given  in  Reference  1 1 . 

A  popular  representation  of  the  regular  radial  function  in  terms  of  spherical  Bessel 
Functions  is  of  the  following  form  described  in  References  9  through  12. 


jeml(c,£)  = 


1  r(£2-01/2 

K,  t 


-r  X  (+i)k+m~‘ 


(2  m  +  k)\ 
k\ 


d:‘jm+k(cZ) 


(22) 


The  above  series  expansion  is  absolutely  convergent.  The  corresponding  series  expansion  for  the 
irregular  radial  function  is  obtained  by  replacing  the  regular  spherical  Bessel  Function  with  the 
corresponding  irregular  spherical  Bessel  Function.  However,  in  the  case  of  the  irregular  radial 
function  the  series  expansion  is  no  longer  absolutely  convergent,  and  it  is  at  best  an  asymptotic 
expansion. 
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The  author  has  studied  the  convergence  properties  of  about  a  half  dozen  different 
representations  of  the  irregular  radial  functions  based  on  integral  representations  of  the  following 
form: 

ne.,(c,4)  =  \dr,K^,mS"](c,r,)  (23) 

Flammer 9  provides  a  general  discussion  for  the  construction  of  integral  representations  of 
the  above  variety.  None  of  these  representations  has  satisfactory  behavior  at  high  frequencies. 
The  high  frequency  convergence  property  of  the  irregular  radial  functions  and  the  angular 
functions  of  the  second  kind  is  due  to  the  fact  that  the  expansion  coefficients  df(c)  represent  an 
alternating  series  of  large  terms  leading  to  numerical  round  off  due  to  the  finite  precision  of  the 
calculations. 
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IV.  SPHEROIDAL  T-MATRIX 


In  the  case  of  the  spherical  T-matrix  description,  the  expansion  of  the  scalar  Green’s 
Function  and  Green’s  Dyadic  in  terms  of  spherical  wave  functions  is  often  used  in  the  derivation 
of  the  spherical  T-matrix  description  from  the  scalar  and  vector  Helmholtz  Equation.  The  vector 
Helmholtz  Equation  in  spheroidal  coordinates  is  non-separable,  and  there  exists  no  closed  form 
solution  for  the  corresponding  Green’s  Dyadic  in  terms  of  spheroidal  wave  functions.  To 
overcome  this  difficulty,  Dr  Hackman  made  use  of  Betti’s  Identity,  which  is  a  generalization  of 
Green’s  Theorem,  to  generate  a  system  of  linear  equations  for  the  scattered  field  in  terms  of  a 
collection  of  basis  functions  of  the  scalar  and  vector  Helmholtz  Equation. 

The  author  introduces  the  following  basis  for  the  regular  and  outgoing  basis  functions  for 
the  scalar  Helmholtz  Equation  in  terms  of  spheroidal  functions: 


^¥am,=  jeml{c,^)Scml{c,Ti,(p) 


Voml=heml(c,£)Saml(c,T],<p) 


(24) 


The  functions  Saml{c,r],q> )  are  the  following  basis  functions  of  the  scalar  spheroidal  harmonics: 


e{m)  (l)  Jcos(w^),£T  =  0 
2xhmI{c)  ml  ’  [sin(«J4?),<7  =  l 


£(m)  = 


j  1  ,m  =  0 
[2,m  *  0 


(25) 


Define  the  following  basis  functions  of  the  vector  Helmholtz  Equation: 


l,  <7  ml 


Vx(a^L) 


V  =— VxF 

V2  ,oml  V^K!  ,<xm/ 

cT 

=  -vrL 

CL 


(26) 
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Here,  Xml  (c)  is  the  eigenvalue  of  the  spheroidal  wave  equation,  and  cT  =  kTf  and  cL  =  kLf  are 
the  dimensionless  wavenumber  for  the  transverse  and  longitudinal  degrees  of  freedom.  The 
functions 


Vo  ml  ~ 

vL^hem,(cL^)saml(cL,T],(p) 


are  the  outgoing  basis  functions  for  the  scalar  Helmholtz  Equation  for  the  transverse  and 
longitudinal  modes.  The  vector  a  is  the  conformal  killing  vector  for  dilatations  whose 
components  are  as  follows: 


a  =  (x,y,z). 


(28) 


This  vector  satisfies  the  following  condition: 

V  a  =  p 

H  v 

The  curl  in  three  dimensions  is  defined  as  follows: 


in  terms  of  the  Levi-Civita  Tensor 


(29) 


(30) 


aaPr  — _ _ _ 

and  the  covariant  derivative. 


(31) 


Define  the  vector  spheroidal  harmonics  for  the  elastic  solid  as  follows: 

Aomt  =  -f=  (Wi-^  -  *7  )s*m,  cp) 

A.oml  =  TX — (^V^- ^  dn+<P  j - -  ^<p)Saml{CT,TJ,(p) 

V4/  V1-7?' 

A.Ornl  =  i^oml  (Ci 


(32) 


In  the  limit  /  — >  0  the  vector  spheroidal  harmonics  approach  the  vector  spherical 
harmonics.  They  obey  the  following  orthogonality  conditions: 


10 


CSS/TR-02/18 


<^d<pdr}Al  (Tml  •  =  Sf8"'Q.Tr 

fadVA2'Cml.A2'a,m,r=d;'S:'Wr 
(§d(pdr]A3aml  •  A3o..mT  =  S°  S™S', 
<^d<pdj]AUanil  •  A3a.mT  =  0 
($d(pdT]A2<yml*AXa,m,r=0 
(jd<pdrjAltTm/*A2a.mr=0 


(33) 


a"r=s'~~rrA . a . .  jdm2sml(cT,n)smr(cT,j]) 

The  vector  basis  functions  VTam,  have  the  following  asymptotic  behavior  in  the  limit  £  — >  +°o ; 

,<7ml 

(34) 

cr 

— *  dgheml(cL,£)A2  , 

CL 

The  above  asymptotic  behavior  of  the  vector  basis  functions  coupled  with  Betti’s  Identity 
leads  to  the  following  important  orthogonality  condition  for  the  vector  basis  functions.  Let  S  be  a 
smooth,  closed  surface,  S „  be  the  sphere  at  infinity,  Fbe  the  volume  whose  boundaries  consists 

of  the  union  of  these  two  surfaces,  and  ( u,v )  be  a  pair  of  regular  solutions  of  the  vector  Helmholtz 
Equation  in  this  volume.  Then  Betti’s  Identity  implies  the  following  surface  integrals  over  the 
closed  surface  S  and  the  sphere  at  infinity  are  equal: 

||<i4{r(w)*v-w*/(v)}  =  || dA{t(u)*v-u»t(v )}  (35) 

^  s. 

By  replacing  the  functions  u  and  v  with  the  vector  basis  functions  and  using  the 
asymptotic  expansions  of  the  vector  basis  functions  and  Equation  33  to  evaluate  the  surface 
integral  over  the  sphere  at  infinity,  we  obtain  the  following  orthogonality  condition  for  the  vector 
basis  functions  on  the  closed  surface  S : 
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ljMmeV,'J.RcV^m,.-f.tV,M  ./(ReF,,..,.))  =0 

5 

J >  •  n-.-.r  -  ^  •  <(R,v-r))  -  o 

s 

\\dA  {/(Re  VT  gm, )  •  VT.  a.m,r  -  Re  Vt  gm,  •  t(Vt,a.m.r ) }  = 


-ifi/krH)  8°a  s: Q“ r  =  r'  *  3 
-/(A + 2//)  /  *t  (-o'-rs:s:s;\z = r  ■ = 3 

0,T5*  r' 


(36) 


The  above  orthogonality  condition  is  the  basis  for  the  spheroidal  T-matrix  approach,  which  uses 
the  above  orthogonality  condition  to  re-express  the  boundary  conditions  of  the  incident, 
scattered,  and  interior  fields  into  a  system  of  linear  equations  based  upon  expansions  of  the 
incident,  scattered  and  surface  fields  in  terms  of  these  vector  basis  functions. 

Consider  the  case  of  acoustic  scattering  from  a  solid  elastic  target  immersed  in  water. 

The  boundary  conditions  between  the  exterior  and  interior  fields  are  given  by  the  following 
conditions  of  continuity  of  normal  displacement,  continuity  of  the  normal  traction,  and  vanishing 
of  the  tangential  interior  traction. 


ii+  •  n  =  u_  •  n 

t+*n  =  t_*n  (37) 

t_y.fi  =  0 

The  single  vector: 

v:m,=-^he„(c4)Sml(.c.ri,<P)  (38) 

c 

is  the  vector  basis  function  for  the  fluid,  which  may  be  regarded  as  the  limit  /i  — >  0  of  the  vector 
basis  for  an  elastic  medium.  The  following  vector  is  the  corresponding  traction  vector  of  this 
basis: 


?(0=^«(v*0  = 

2wk 2  (39) 

- - hem,{c4)Saml{c,v,(p)n 

c 


Make  the  following  expansions  of  the  incident  and  scattered  field  in  the  fluid  in  terms  of 
this  basis: 
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«»,=2>.-ReC 

ami 

"s-^L/.jC  (40) 

ami 

Applying  Betti’s  Identity  to  the  total  field: 

UT</laI  =  Uinc  U scat!  ~  ^  a<yml  ml  0 

ami 

and  the  regular  and  irregular  basis  functions  on  the  surface  S  of  the  scatterer,  one  arrives  at  the 
following  linear  equations  between  the  expansion  coefficients  of  the  incident  and  scattered  fields 
and  the  surface  fields 


=-'-4  \m,  *C  *'(Oi 

^  S 

=  +‘jw  •  ReC  *'(ReOi 


(42) 


Using  the  fact  that  traction  in  the  fluid  is  proportional  to  the  normal  vector  to  the  surface 
and  the  continuity  of  the  normal  displacement,  one  can  rewrite  these  equations  in  the  following 
form: 


=-'4r  J 

*  ‘  (43) 

/-  =  +/-£-  • »)« •  Re  C,  -(«  • »)« •  dRe  C)} 

Next,  one  may  apply  Betti’s  Identity  to  the  regular  vector  basis  functions  of  the  solid  and 
the  surface  displacement  and  traction  of  the  interior  fields  to  arrive  at  the  following  equations. 

\dA{t_  .ReK, ./(ReR,  „,)}  =0  (44) 

S 

Use  the  boundary  conditions  at  the  surface  to  replace  this  integral  with  the  following: 

•»)(«•  Re •'(ReR,itr„))  =  0  (45) 

S 

Now  make  the  following  expansions  in  terms  of  the  scalar  harmonics  in  the  exterior  and  the 
vector  harmonics  in  the  interior: 
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t+*n  =  ^0!amlSaml  (c,  rj,  <p) 

ami 


«-  =  Z  Pr#mlA r. 


ami 


rami 


(46) 


Substituting  the  above  expansions  into  Equations  43  and  45  one  arrives  at  the  following  linear 
system  of  equations: 

a  =  -(Q/5-Ma) 


f  =  +  ReQfi-ReMa 


0  =  Pa-Rp 


Q, 


A, rV 


=  -i -jjw  \dA{{n •  t(V* )(n •  Ar  J} 


n,n 


jdA{(n»V")Sn.) 


(47) 

R,„r,=  jdAmeV,,)'^,} 


The  indices  crml  have  been  grouped  into  the  single  index  n  in  the  above  equations.  The 
solution  of  the  above  systems  of  equations  is  given  by  the  following  T-matrix,  which  relates  the 
incident  and  scattered  fields  in  the  fluid. 

T  =  -(ReQR-'P-ReM)(QR-'P-MY'  (48) 

The  above  T-matrix  is  Method  B  described  in  Reference  3.  This  is  the  method  used  to  perform 
the  calculations  in  the  next  section. 
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V.  SAMPLE  CALCULATIONS 


Figures  1  and  2  depict  a  comparison  of  the  results  of  the  present  author  and  Dr.  Hackman 
for  the  scattering  from  a  4: 1  solid  aluminum  prolate  spheroid  at  end-on  incidence  (0  degrees)  and 
30  degrees  incidence. 
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Comparison  of  Scattering  from  4:1  Aluminum  Spheroid  at  30  Degrees 


1  3  5  7  9 


KL72 


FIGURE  2.  FORM  FUNCTION  AT  30  DEGREES  INCIDENCE 

Calculations  made  by  the  present  author  and  Dr.  Hackman  agree  within  machine 
precision  in  the  above  two  cases. 
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Figure  3  depicts  a  polar  plot  of  the  form  function  at  kL/2=5.5,  which  corresponds  to  the 
frequency  of  the  large  flexural  wave  depicted  in  Figure  2.  The  directivity  pattern  of  the  flexural 
resonance  at  kL/2=5.5  has  the  characteristic  shape  of  a  quadrupole. 


FIGURE  3.  FORM  FUNCTION  AT  KL/2  =  5.5 
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Figure  4  depicts  a  polar  plot  of  the  bistatic  form  function  at  kL/2=5.5  at  end-on  incidence. 
The  forward  scattering  is  strongly  peaked  at  the  radiation  lobes  of  the  flexural  wave. 


FIGURE  4.  BISTATIC  FORM  FUNCTION  AT  KL/2  =  5.5  FOR  END-ON  INCIDENCE 


21/22 


CSS/TR-02/18 


Figure  5  illustrates  a  two-dimensional  rendering  of  the  form  function  of  a  4: 1  aluminum 
spheroid  as  a  function  of  frequency  and  aspect  angle.  The  first  three  flexural  resonances  of  the 
spheroid  are  clearly  legible. 


L/D=4  Aluminum  Prolate  Spheroid 


kL/2 

Aug  16  13:24 


10.00 


-40.00 


FIGURE  5.  TWO-DIMENSIONAL  PLOT  OF  THE  FORM  FUNCTION  AS  A  FUNCTION  OF 

FREQUENCY  VERSUS  ASPECT  ANGLE 
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APPENDIX  A  -  COMPUTATION  OF  EIGENVALUES  AND  EIGENVECTORS 


This  appendix  describes  the  computation  of  the  eigenvalues  and  coefficients  df  of  the 
spheroidal  wave  equation. 

The  expansion  coefficients  df  satisfy  the  following  recurrence  relation: 

«;< 2 + (  p:  -  k,  k‘ + = o 

„  _  (2m  +  k  +  2){2m  +  k  +  l)c2 
k  (2m  +  2k+3)(2m  +  2k  +  5) 


=(m  +  k)(m  +  k  + 1) 

2  2(m  +  fc)(m  +  £  +  !)-  2m2  - 1 
+C  (2m+2k-l)(2m  +  2k  +  3) 


rk  = 


k(k-l)c2 


(2m  +  2k-  3)(2m  +  2k  - 1) 


Define  the  ratio 


iml 

N  =-am  a-  — 

1Sk  -  Uk- 2  jml 
“k-2 


These  ratios  obey  the  following  forward  and  backward  recurrence  relations: 

<07 


N, 


N  _ 


^2+6  ~  fis  A, 


S  =  (l-m)  mod2 
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The  eigenvalue  A  is  chosen  such  that  the  following  limit  vanishes: 
lim,^  Nk=0 

The  equation  for  the  eigenvalues  is  obtained  by  using  the  forward  recurrence  relation  to 
compute  t  he  value  of  the  above  ratio  at  k  =  2+1  —  m .  Next,  one  chooses  a  sufficiently  large 
value  of  k  =  K  such  that  the  ratio  NK+2  =  0  vanishes  and  computes  the  value  of  the  ratio  of 
k  =  2  +  1  -m  from  the  backward  recurrence  relation.  The  eigenvalue  is  adjusted  to  make  the 
difference  between  these  two  estimates  of  the  ratio  at  k  =  2+1  —  m  vanish. 

Once  the  eigenvalue  has  been  obtained,  the  values  of  the  above  ratio  for  k  <2  +  1 -m  axe 
computed  from  the  forward  recurrence  relation.  For  indices  k>2  +  l-m,  the  above  ratio  is 
calculated  from  the  backward  recurrence  relation.  The  rationale  for  splitting  the  above 
calculation  into  a  combination  of  forward  and  backward  recurrence  is  to  avoid  loss  of  precision 
by  recurring  in  the  direction  of  increasing  magnitude  of  the  ratio.  The  expansion  coefficients 
df  are  obtained  from  the  recurrence  relation: 


dms'=  1 


jml  _™ 
“  k+2  ~~ 


N, 


k+2  jml 


cc 


The  above  prescription  defines  the  expansion  coefficients  df  up  to  an  arbitrary  normalization 
condition.  One  adopts  the  normalization  condition  of  References  1  and  2  defined  below: 

5^/(0)  =  /f  (0),l-m  even 


S™  '(0)  =  P,m  '(0),  1-m  odd 

X',_Un  (2w  +  2m  +  2£)!  /  -  y  i-m-S)/2 _ (l  +  m  +  S)\ _ 

2J<  l>  2 2n+\n)\(n  +  m  +  8)\  2n+s  2^m \(l  +  m  +  d.) i 


8  =  (/-w)mod2 

An  alternative  method  of  computing  the  eigenvalues  and  eigenvectors  of  the  spheroidal 
wave  equations  utilizes  the  fact  that  the  recurrence  relation  for  the  expansion  coefficients  can  be 
represented  in  the  following  matrix  form  as  an  eigenvalue  problem: 
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'A 

as 

0 

o  > 

>m/  ^ 
“6 

f  dml  ^ 
a6 

Ps+2 

0 

jml 

“8+2 

=x 

jml 

“8+2 

0 

ys+a 

ft 6+4 

«<S+4 

dml 

US+4 

jml 

“8+4 

,  0 

0 

y 6+6 

A+6> 

jml 

\as+ b) 

iml 

\  “8+6  J 

The  matrix  can  be  transformed  into  a  symmetric,  tri-diagonal  matrix  by  defining  the 
following  expansion  coefficients: 

n0  a  1,  n  =  0 


j-j  _  ^£+2(n-l)  j-j 


7s+2n 


n- 1 


(X; 


'6+2  j 


7 5+2(j+\ ) 


n>  0 


en  -df,  n  =  0 


en=nX,'+*>  n  >  0 

These  coefficients  obey  the  recurrence  relation: 


an  —  \j^6+2nYsH 


S+2n+2 


S+2n 


a0e]  +(b0-A)e0=  0,n  =  0 


anen+i+(b„-Z)e„+an_lbn_ , 


=  0,n>  0 


The  eigenvalues  and  the  eigenvectors  are  solutions  of  the  eigenvalue  problem  for  the 
following  symmetric  tri-diagonal  matrix: 


ib° 

ao 

0 

0  N 

^0 

eo 

ao 

b\ 

0. 

e\ 

=  A 

0 

a. 

«2 

*2 

e2 

0 

«2 

*3  J 

ve3y 

^e3> 
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This  method  of  solving  for  the  eigenvectors  is  numerically  more  robust  than  solving  the 
transcendental  equation  for  the  eigenvalues  presented  previously.  In  addition,  this  method  is 
readily  generalized  to  the  case  the  wavenumber  (c)  is  complex.  Note:  the  user  is  required  to 
solve  the  above  eigenvalue  problem  for  even  and  odd  parity  eigenvalues  separately. 

The  value  of  the  expansion  coefficients  df  for  —2m  +  5  <k<— 2  +  5  can  be  calculated 
from  the  backward  recurrence  relation  for  the  ratio  of  adjacent  coefficients. 

<-lS— 

*  M-*~Nk+2 

5  =  (l-m)  mod2 

Alternatively  one  can  use  the  following  forward  recurrence  relation  for  the  coefficients 

d_k  • 


'1,  A  +  «r.M) 


r-k 


iml 

— y—,k  <  2m  +  2-S 

iml  7 
<  “~k+2 

0 ,k  >  2m  +  2  —  S 


Since  the  term  r_2m_2+s  =  0  vanishes,  the  above  recurrence  relation  expresses  the  ratio  as  a  finite 
continued  fraction.  The  expansion  coefficients  dZ[  are  given  by  the  following  expression: 

(n;,o  r_2n_2.  +s)df  =  dml_2+s,n  <  m 

In  the  case  k  >  -2m  -  2  +  8 ,  the  coefficients  d'fp  have  a  simple  zero  in  the  limit  p  ->  0 . 
In  this  case,  one  is  interested  in  the  coefficients  d^p  for  k  >  -2m  -2  +  5  defined  by  the  relation. 

iml 

One  may  use  the  following  forward  recurrence  relation  for  the  coefficients  d™‘p : 
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f.„  /(«  -A+JC'-.,.,),*  >  2m  +  2-5 


iml 

k  =  2m  +  2-S 

iml  ? 
a-2m+£ 

jml 

— ^-,&  >  2m +  2 -S 

iml  7 
.  “*-2|p 


a”* ,k  *2m  +  2-S 


(2w-l)(2m  +  l) 
c2 

(2m-3)(2m-l) 
Here  the  coefficients  ^ 


,k  =  2m+2— 8,8  =  0 
A:  =  2m  +  2-S,S  =  1 

are  given  by  the  following  expression  in  terms  of  the  ratios 


01”=  0  2(fw+w)— 2+<5  'l^—2m+S 


jml 

a+2  (m+n)+2-8\p 
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APPENDIX  B  -  EIGENVECTOR  ORTHOGONALITY  RELATIONS 


This  appendix  contains  the  orthogonality  relations  satisfied  by  the  expansion  coefficients 

iml 

dk  . 

The  expansion  coefficients  df  satisfy  the  following  relationships: 


y  2 

(2m  +  *)!  ,  sr 

tSi(1m+2k+\) 

ak  ak  - 
k\ 

£«/a„,= 

l-m 

(2wi  4~  2k  + 1)  k !  £k 

2  (2 m+k)\  k 

A.,=  P'7S2fe'7)«’fe’?) 

-1 

The  above  relationships  represent  the  orthogonality  and  completeness  of  the  expansion 
coefficients  viewed  as  eigenvectors  of  an  infinite  dimensional  tri-diagonal  matrix. 
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APPENDIX  C  -  GREEN’S  FUNCTION  EXPANSION 


This  appendix  contains  the  expansion  of  the  scalar  Green’s  Function  and  the  plane  wave 
in  terms  of  spheroidal  wave  functions: 


G(r,r') 


exp[+/&  |  r  -  r '  |] 
An\r-r'\ 


=  »*£  heml  (c,  &  )jeml  (c,  ^)Saml  (tj,  (p)Scm,  (77  \  <p  ’) 

ami 


exp[+/£  •  r  ]  =  exp[+/Ar  cos(/)]  =  4;r]T  (+*')'  jem,  (c,  %)Saml  (*7>  (P)sami  07 

ami 


cos(7)  =  cos(t?)  cos(t?') + sin(t?)  sin(t?')  cos  (<p  -q>') 
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APPENDIX  D  -  TRANSFORMATION  FROM  SPHERICAL  TO  SPHEROIDAL 


This  appendix  contains  the  transformation  between  spherical  and  spheroidal  wave 
functions. 

The  transformation  between  spherical  and  spheroidal  wave  functions  is  given  by  the 
following  transformation,  where  the  matrix  Baml  a.m.r  is  a  unitary  matrix: 

^ml^aml  ~~  ^  ^ ami ,a' o' m'V 
G'm'V 


Wornl  ~  B 


G'm'V 


G'm'V 


R  y  (+0  ‘‘1  2  (r+w)!  , 

ami,a'mT  a  m  ]j  (2/'+l)  (/>- m)! 


d-i  _  (+f')/  [  2  (/  +  /w)! 

ami.awr  a  m  \  (2/ + 1)  (/ -  ill) !  ^ 


_  Dt 

^  GmltG'm*V  ^  Gml,G'm'V 
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APPENDIX  E  -  COMPUTATION  OF  RADIAL  FUNCTIONS 

This  appendix  describes  various  representations  of  the  radial  functions.  These 
representations  are  based  on  integral  representations  of  the  following  form  described  by  Flammer 
in  Reference  1 . 

*Jc,£)  =  )w(lvK\c,ri) 

a 

The  following  Theorem  and  Lemma  are  from  Reference  1 : 

THEOREM: 

Let  Z,£  and  Ln  be  the  following  differential  operators: 

Let  be  a  solution  of  the  following  partial  differential  equation  in  the  complex  domain  D: 

(L§-Ln)K(t,i»  =  0 

Let  [a,b]  be  the  contour  in  the  complex  domain  D  such  that  the  following  integral  vanishes: 

T1  )-m  TJ)}  = 

a 

a  -  n1 )  {0„K(f .  v»Oc,  n)  -  (c,ml  =  o 

Then  the  following  contour  integral  is  a  solution  of  the  radial  equation  in  the  complex  domain  D : 

Rm,(c,Z)=)dt,K(Z,n)s!!!(c,i}) 

a 

(Lg  +  Ami  )Rm,  =  0 
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LEMMA: 

Suppose: 

HZ’  H,  (p)  =  K(Z,  TJ)  exp  [+inup] 
is  a  solution  of  the  Helmholtz  Wave  Equation. 

(V2  +k1)y/  =  0 

Then  the  function  K(^,  rf)  is  a  solution  of  the  partial  differential  equation: 

(LrL„)K(t;,r1)  =  0 

Some  of  the  examples  of  kernels  satisfying  the  above  partial  differential  equation  listed  in 
Reference  1  are  listed  below: 

(f-ir/2(l-772r/2exp[+/c^] 

ym(CV(<f-i)(i-772)) 

ZvJm(cyl(?-m-ri2)) 

Jm  (cy[(£?  — 1)(1 — Z72 )  sin(or))  exp  [tic^Tj  cos(a)] 

jm(c{Z±me-^)ml\^2r\c^±n)r 

h'\ (C(t;  ±  mz2  -  i)ra/2  (i  -  n2  r1  «z  ±  mr 

The  classic  expansion  '"2  of  the  radial  functions  in  terms  of  spherical  Bessel  Functions  follows: 


E-4 


CSS/TR-02/18 


•  l  g^  1  (£ 2  ~  0"'  ~  V  /  .  --,2 n+S+m-l  (2m  +  2 «  +  <?)!  ,m/  /  C\ 

Jeml(C>?)  ~  j^r)  gm  JLj  ^  '  (2n  +  S)\  2n+sJ2n+m+S\  £) 


/  _  g-.  j_(£  j)  'ST' ,,j\2n+S+m-l  (2m+.  2«  +  <?)!  ■„,/  /_£\ 

^m/W»§)  ^(r)  '  (2n  +  S)\  2”+<*  2n+m+^  ^ 


8  =  (/- w)mod2 


Afr>  y(2«  +  2/»  +  <?)! 


jml 

“ln+8 


The  expansion  for  the  regular  function  jeml  is  absolutely  convergent,  whereas  the 
expansion  for  the  irregular  function  nem,  is  an  asymptotic  expansion  unsuitable  for  radial 
coordinates  near  unity. 

The  expansions  of  the  radial  functions  in  terms  of  the  analytic  continuation  of  the  angular 
functions  are  given  by  the  following  expressions: 

Kml 

Kml 


^(2m  +  2n  +  8)\  ^ml 


(2m  +  2S  +  l)(l  +  m  +  S)\ 

+  m  +  (2  n  +  8)\ 


2  n+8 


«S’=+- 


2'~(2m)'(— )!(— )W" 


2m  (2m  +  2n)!  ,m/  ..  . 

- - - - - ; - )  - - d2ll,  even(l-m) 

(2m  -  l)m !(/  +  m)  Ic”-  %&  (2  «)! 


2l~m (2m) !(-  -  •) !(/  +  m  +  1) ! «/? 

r<2)  _ _ 2 _ 2 


2m+i  ~  (2m  +  2n  +  <5)! 


'■ml 


(2m  -  l)(2m  -  3  )m !(/  +  m  + 1)  !c‘ 


m-2 


(2n  +  8)\ 


dZs>  °dd  (1-m) 


The  above  representation  of  the  irregular  radial  function  nem,  in  terms  of  associated 

Legendre  polynomials  is  useful  for  calculating  the  radial  function  for  arguments  in  the  range 
+1  <  £  <  1 .5 .  In  the  above  expressions,  one  uses  the  analytic  continuation  of  the  associated 
Legendre  polynomial  to  the  complex  plane: 
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pr<,z)=(z2-i  rn(4-rm 

az 

e-(r)=(rj-ir'jA"ew 

az 


The  following  expansion  of  the  regular  radial  function  in  terms  of  spherical  Bessel 
Functions  is  sometimes  useful  for  arguments  near  unity: 

je„  =^Z(+i)"“'<C,(0)y„»(cJfM)  (even  1-m) 

X,  =  .„,L  ,i  ,  (odd  1-m) 

°mlly)y]s  - 1  k 

There  is  an  error  in  Reference  1  in  the  case  l-m  is  odd  in  the  above  expansion. 

The  radial  functions  have  the  following  expansions  in  terms  of  the  product  of  spherical 
Bessel  Functions: 

£=cosh(///2) 

Even  1-m 


he„  =  (f: 1  - 1)”'!  t  (2n + 2»  +  1)B («*"2  "2 ) 

rt=0 

,  (+/)'-n,2w(n  +  2/n)!y.  ,  (-«)t(n  +  2m  +  l)t(m+l/2)t  r(m  +  3/2) 

"  d"‘cmm\n\  h  21  k\  Tim +  3 1 2  + 2k) 

x3F2  (— n  +  k,  n  +  k  +  2m  +  \,k  +  \/2;m  +  k  +  \,m  +  2k  +  3/2:+\) 

(+/y-m  (2m  +  l)2m~l  ( n  +  2 m) !  ^  ,  (-ft),  (ft  +  2m  + 1),  (m  + 1  /  2),  T(m  +  ft  -  A:  + 1  /  2) 

d"'cmm\n\  h>  2k  k\  T(m  +  n  +  k+3l2) 


E-6 


CSS/TR-02/18 


Odd  1-m 


heml  =£(?  - l)m/2 £ (2«  +  2m  +  3)^„+1  (^+^/2 )y_„+1  (C^/2 )iC»r 

n=0 


m  (+Q1— 1  (2m + 3)2*(2m  +  n + 2) ! ^  i  (-n),(»H-3/2),(«  +  2m  +  3),  r(m+3/2) 

"  d?cm*'m\n\  h,  *!  Hm+5/2  +  2*) 


'xiF2(-n+k,k +3 /2,n  +  2m  +  k+3\m+k  +  2,m  +  2k +5 /2:+l) 


The  Pochammer  and  generalized  Hypergeometric  function  are  defined  below: 

(*),=r(x+t)/rw=n(x+n) 

n=0 


3  F2(a],a2,a3;bl,b2:z)  = 


y1  (a\)n(a2)n(ai)n  = 

n\(bx\ib2)n 


REFERENCES: 

1.  Flammer,  C.,  (1957),  Spheroidal  Wave  Functions ,  Stanford  University  Press,  Stanford  CA. 

2.  Hodge,  D.B.,  (1970),  Eigenvalues  and  Eigenfunctions  of  the  Spheroidal  Wave  Equation, 
Journal  of  Mathematics  and  Physics,  V.  1 1  pp.  2308-2312 
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APPENDIX  F  -  CONNECTION  IN  SPHEROIDAL  COORDINATES 


This  appendix  contains  a  tabulation  of  the  non-zero  terms  of  the  linear  connection,  its 
derivatives,  and  the  function  P^y  =  daT%,  +T^aT^y  in  prolate  spheroidal  coordinates.  These 

functions  are  used  to  compute  the  second  and  third  order  covariant  derivatives  of  the  scalar  basis 
functions  used  in  the  construction  of  the  vector  basis  and  the  stress  tensor: 


rc  gq  zltl 

*  (£2-l)(f-772) 

_  77(£2-l) 

m  (l-772)(^-772) 


T) 


rt  =  r£  = _ 

*  *  (f-//2) 

_  ■p'?  —  i 

*  to  (f-772) 


r*  _  flf-P 
w  O-^Xf-tf2) 

r,  _  ^q-^2) 

*  (f-l)(f-772) 


= _ jZ 

"  w  (1-772) 

pp  _  pp  —  I  ^ 

*  ft»  +(f-l) 


pf  ^2-l)(l-772) 

w  (f-772) 

r„  _  ,  77(l~?72)(f -1) 
w  (f-^72) 


q-^2)  {1  2? 


2? 


^  p£  _ _ _ 

*  *  tf2-lX^2-*72jl‘  (l2-D  (?~V2) 


} 


2^ 


3  H  — _ 

"to 
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VI,  =3^=+^^- 


a  r«  =9  r*  =_i£±i> 

Vfr  V»{  ^2_^2-)2 


*)  r{  _,  1  fi  ,g2  ,  2£2(£2-l) 

+  „  V  2\  ^  y  f*2  __2\  > 


*  w  (1  -V2)(?-T12) 


(F-tf) 


S  2^(f-l)  1  1 

7  w  (l-^2)(f-72)Xl-/72)  (f-772) 


-V  ,  (l-?72)  ft  oK2  2£2(£2-l), 

V^=+—  —  fl-3g  +— - — } 


(F-v ) 


(F-n2) 


a  ^ 


a  rf  -a  r'  -  +  0 

°5‘  rC  ar*  (  =2  _j.; 


V*=V*.=0 


VJ,  =3{rj,  =o 


^  p?>  _  ^  p?  o+yy ) 

*V  W  CV  W  ^_^2 
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2^/y(l-772)  1  1 

e  *  (f-l)(f-77V(!2-l)  (f-/72) 


a  r„  1  n,„ 2,  27fQ-7f) 

s,fl  (f-ixf-’l2)1  n  (f-V) 

2^/7(l-772)2 
*  w  (<f-?72)2 


aJT..+j£z«L{1_a,>+^ii2, 


V  <P<P 


(g-rt) 


ta  ,  3gflz2!) 

W  (f-l)2(f-^72) 


/«  =+ _ & _ 

**  (f-lX^-772) 


pS  3^(1  ~^2) 

**  (f-l)2(f-^2) 


pt  _  p#  _  _| _ & _ 

w  (f -l)(f -772) 


p>7  =  p7?  = _ 2 _ 

**  (l-/72)(^2 -?72) 


pf  =PS  1 _ 

^  wf  (f-IXf-T/2) 


^7 


O’?  _  D7  _  . _ 

(l-^2Xf-772) 


£2 


n  3&X?-1) 

*"  (1-'7,)j(#2-'7!) 
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pt  =+ — in - 

nm  (1-V2)(t2 -v2) 

3772(f-l) 

nm  (1  -T]2)\¥-712) 
JA  =_£0z rf) 

(?-j f) 

ps  _ 1  Mzff!) 
w  (f-^2) 


.  MlzB 

^  (f-772) 


w  (£2  -  772) 


pi  —  pi  _  _ 
«0?¥>  w? 

=  - 

m  wn 


pn  —  pn  _ _ 


1 

(f-D2 

& 

(l-772)(f-l) 

to 

(l-/72)(f-l) 


pn  —  pn 

<pn<t>  m 


1 

(l-^2)2 


**  =  **  _  £Ozg!> 

££?  ft»?  (g2  —rf) 

pv  =p<p  _  ,  -1) 

fw  #«”7  _^2) 

n*  =  p*  - 1  Mz£> 

ni9  nrt 


p<p  _  p<p 

nm  n<pn 


/72(f-0 

(f-72) 
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p<p 

rv4n 


_  p<p  _  _ 

“  rwff  ” 


p<p 

m 


1 

(i-n2 


(l-772)(f-l) 


)2 


F-7/F-8 


CSS/TR-02/18 


APPENDIX  G 

VECTOR  BASIS  FUNCTIONS 


Approved  for  public  release;  distribution  is  unlimited 


G-l/G-2 


CSS/TR-02/18 


APPENDIX  G  -  VECTOR  BASIS  FUNCTIONS 


This  appendix  contains  the  components  of  the  vector  basis  functions  and  their  covariant 
derivatives  for  an  elastic  medium  in  spheroidal  coordinates.  One  may  use  the  index  n  to  signify 
the  group  of  indices  ami .  Here,  y/Tn  and  y/Ln  are  the  scalar  basis  functions  for  the  transverse  and 
longitudinal  modes.  The  vector  basis  functions  are  defined  by  the  following  relationships. 

<*n =/ 2v 

^=0 


g  =  Det(gftv)  =  f6(?-n2)2 


yJS 


+2  vy  +aavavy„ 


fry  __  t  S^Sva^^ pWl 

'  Jg 

+k2ravvyn  +Wfjvy  +  aaVaVMVyn 

V„(F,>),  =7-V,V„tf 

kL 


First,  one  tabulates  the  components  of  the  second  and  third  order  covariant  derivatives  of 
the  scalar  basis  functions.  These  functions  are  used  to  construct  the  vector  basis  and  its  covariant 
derivative.  Since  the  Riemann  Curvature  Tensor  vanishes,  the  covariant  derivatives  commute  in 
a  coordinate  basis.  Thus,  the  higher  order  covariant  derivatives  of  the  scalar  basis  functions  are 
symmetric  tensors.  The  radial,  angular,  and  azimuthal  wave  equations  have  been  used  to 
eliminate  second  and  higher  order  derivatives  with  respect  to  a  given  variable: 
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W-rf)  4  n  0-772)  ^  1  m  JiJ,  w 


v^-+o=S^^+^-^p‘r-^-cV-o$jIr 


'W = a  Ar  +(pi^yfofyy "  ^ 


v«v#r=aA,r-^r^d,r 


v,v#nr»a^+ji^a,r 


VyV>y  =  -«V+(1  (^-^)  1}[^-^] 


V.V.V^r  =  -  e23;7(1  Tp-  2 

5  ^  (|2-l)(f-772)  5  ^ 

2^7y(l-772)  2  1 _ (I-772)  3 

+  (^2-l)(f-^2)L(f-772)  +  (f-l)  (f -l)(^2-^2)j  ,r 


3(l-772) 

(f-l)2(f-/72)2 


[^2+^+^V-3^4]3^ 


+— ^ — [A + 6-c2^2 + +— ja^ 
(f-D  (f-i)  ^ 


4  rl  _2e2i  .  6/«2  3{\.  —  Tj2) 


{2c2+-f—[A-c2f2]+ 


(?-»'  (f-1) 


Of -l)2  (f-1  )(f-/72) 


-[W<f  + 


(f-0 
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.  2{ri(f!-l)  ,  2  ,  1 

T  -i  ^  |_  -j  > 


(f-D 


O-^X*2-*2)  (f-72) '  a-^2)  a-n2)(Z2-v2) 

3(f-l) 


1 


d-772) 

V 

d-^72) 


{/L  +  6-c  77  “ 


[^+7j2+^rj2-3^]dv¥ 
(m2  +8) 


d-^72) 


{2c2 


[X-c2if}+ 


(i-r> 


6m 


3(f-l) 

,2\/e2 


-[>1  — C2772  — 


d-^2)2  a-n'xr-i) 


m 


a-*r) 


•]}y 


= +-3(1(j?  j^2)  -  (™2 + 2>9^ 

W>=-£-^—  +  2 


o-^2) 


V 


(f-1)  (f-/?2)  (f-ixf-^2) 

.2*2  _2«2  .  2m2  m2 


(f-iXf-V) 

(f-lXf-/?2)2 


[31-2c2^2-cV  + 


(£-1)  (l-//2) 


[3-4|2  +77 


1  “2 


(f-D 


[i-c2|2  + 


rR.V' 


1 


(f-1  X^-172)2 

2  2 


(f -1)J  ,V 

[3t74  -  3//2  +  ^  V  -  3£2  +  2£4  ]9_^ 


^V,V^  =  77[- 


(|2-1) 


a-??2)  (f-^2)  o-^xf-r) 

£  ,m  ,  2  2  2e2  2m2  m2  . 

{3A-2cV  -c2f  }y, 


(l-tfXZ2 -TJ2) 

(1  -tfyf-rff 


(I-772)  (f-1) 


[3-4t72+^2]3„^ 


1  r1  22  m2 


a-r) 


i 


(1-T72X£2-T72)2 


a-*r) 

[3£4  -3f  +£V  -3^2  +2t74]3^ 
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^v,7v?)^=a^,a^+[ 


(l-*72) +  (? -772)33^  [(£2  -1) +  (f  -WdM 


.r _ & _ + _ §1 _ + _ _ p  y, 

L0-*72Xf-772)  (l-72)(f-D  (f-l)(f-72)  ^ 


= +  _(11}  J2  ^[gAy-»?3Ay] 


{?-!)' 


(f-iy  (f -i) 


(f-D  r 

,,  w2  +  2 


vvv„  7(1-V2)(#!-1)-,  s 
V«V„V^  = - 2_[^) - 

«2->?2)2  ,r  tf2-<72)2 
W-^)„  cy  «* ...w+  ^  „ 

(f-v2)1  4  «2-ir  «2-Dr 


-- 11 


v,v,v,t/=+fi!J^-Saia„r 

|tf-‘)r.i  „y _ g!_w_  ^  w 

V-W  n  Q-i'r  (1 


G-6 


CSS/TR-02/18 


The  component  expression  for  the  above  vectors  in  terms  of  the  scalar  basis  functions  are 
as  follows: 

y.  (K ), = '@y  -  w>y.  i 


Kjkvy =crW,  +*y.  +§ w 


K4Ky\=c2Tw: +2vx +|^^vtv,y: 


*tJ\(k%=2  vx +|mV|vx 


(0;=-3^,‘ 

CL 

y,\=—^y 

CL 

(rx=-zy 

CL 
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The  components  of  the  covariant  derivatives  of  the  vector  basis  functions  are  as  follows: 

:VIVe(K.l){=-H70 -J?2)VjV,r.r 


^r.\  - w  V-ix£f)vx 


j-JX7,<r.),  =  -1)VX  -tfi(~m^p[-/V{7X  -£V„VX  ] 

jVIv,(K.')„  =  +(£2  -1)VX  -^yv»VX 


j  VXv,(k;  ), = -  FSS.  ] 
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kr&V'rt),  +3V{v{rI  +|f^' Wstf 


(f-72) 


(f  ->7  ) 


*r&?,<r.% =<*  +3v„v,tf  t^iw/  +|^  WX 

*rVW?)f =4( f! -Dd-^)^: +3v,vx +|liiliv;v,v,< +|^v,v,v„r; 

*rV4W!),  =4nvy  +3v{v,tf  w,rJ  +|I ^|wx 

kr4\v,<.v;h = +3vsvx +|^v,w  +|^  wx 
kT^e(r.\  =3  v{v,tf  +^|wX 

trJZWX  =  3V„VX  ^V(T,V/  +|^W,^ 

K^,(n\ =4&y  +3v{vx +§r-ipy  wx  *|^wx 


kr  #,K')rW +3V.VX +f£^  WX +§c^  WX 
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VOs  =  — V;Vsr.‘ 

CL 

v„K3),=-v,vX 

Cl 

v  (y3)  =— v  v  wL 

T  p\r  w  /?>  T  p  T  pr  n 

CL 

Vj(0„  =^w.‘ 

CL 

W’)(=-W 

CL 

'WV-VjVf 

CL 

v,(n3){ =— v{vx 

CL 

W)»  =— v,vx 

Cl 

VO,  =^v„v,r‘ 

CL 
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APPENDIX  H  -  PARITY  OF  BASIS  FUNCTIONS 


This  appendix  tabulates  the  parity  of  the  various  functions  and  their  derivatives  with 
respect  to  the  transformation  i]  and  the  dependence  of  the  functions  on  the  first  order 

derivative  with  respect  to  the  azimuthal  angle.  These  symmetries  are  important  for  determining 
the  non-zero  matrix  elements  obtained  from  the  surface  integrals  for  a  symmetric  target  such  as  a 
prolate  spheroid  or  finite  cylinder. 

Table  H-l  contains  the  polar  and  azimuthal  parity  of  the  vector  harmonics.  Table  H-2 
contains  the  parity  of  the  components  of  the  vector  basis  functions,  and  Table  H-3  contains  the 
parity  of  the  stress  tensor  formed  from  the  covariant  derivative  of  the  vector  basis  functions.  The 
interpretation  of  the  polar  parity  column  is  as  follows:  +1  signifies  the  function  transforms  as 
F{-rf)  =  +(-l)'“mF(/7) ,  -1  signifies  the  function  transforms  as  F(-rj)  =  -(-1  )‘~m  F{if) ,  and  0 
signifies  the  function  is  zero.  The  interpretation  of  the  azimuthal  parity  column  is  as  follows:  Y 
signifies  the  function  is  proportional  to  an  odd  number  of  derivatives  of  the  function  with  respect 
to  the  azimuthal  order,  and  N  signifies  it  is  proportional  to  an  even  number  of  derivatives. 

The  significance  of  the  azimuthal  parity  is  that  it  denotes  the  coupling  of  different  parities 
in  the  azimuthal  integration.  The  significance  of  the  polar  parity  term  is  that  it  dictates  the  values 
of  order  /  for  which  the  integrals  over  the  coordinate  -1  <  7)  <  +1  are  non-zero. 

TABLE  H-l.  PARITY  OF  VECTOR  HARMONICS 


VECTOR  HARMONIC 

COMPONENT 

POLAR  PARITY 

AZIMUTHAL  PARITY 

First  Vector 

1 

£ 

0 

0 

1 

+1 

Y 

i 

<p 

-1 

N 

Second  Vector 

2 

0 

0 

2 

n 

-1 

N 

2 

<p 

+1 

Y 

Third  Vector 

3 

S 

+1 

N 

3 

V 

0 

0 

3 

(p 

0 

0 
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TABLE  H-2.  PARITY  OF  VECTOR  BASIS  FUNCTIONS 


1  VECTOR  BASIS  FUNCTION 

COMPONENT 

POLAR  PARITY 

AZIMUTHAL  PARITY 

First  Vector 

i 

4 

-1 

Y 

i 

V 

+  1 

Y 

1 

<p 

-1 

N 

Second  Vector 

2 

4 

+  1 

N 

2 

V 

-1 

N 

2 

9 

+  1 

Y 

Third  Vector 

3 

4 

+1 

N 

3 

V 

-1 

N 

3 

9 

+ 

Y 

TABLE  H-3.  PARITY  OF  STRESS  TENSOR  FORMED  FROM  VECTOR  BASIS  FUNCTIONS 


VECTOR  BASIS 
FUNCTION 

COMPONENT  OF  STRESS 
TENSOR 

POLAR 

PARITY 

AZIMUTHAL 

PARITY 

First  Basis  Function 

4) 

-1 

Y 

(£,77  )(//,£) 

+1 

Y 

{£,<P)(<P,€) 

-1 

N 

-1 

Y 

i 

(n,<p)(<pji /) 

+1 

N 

1 

( <p ,  <p) 

-1 

Y 

Second  Basis  Function 

2 

<f.f) 

+  1 

N 

2 

(£>77)(*7>£) 

-1 

N 

2 

(€,<p)(<p,Z) 

+  1 

Y 

2 

W,V) 

+1 

N 

2 

(TJ,<p)((p,TJ) 

-1 

Y 

2 

(<P,(P) 

+  1 

N 

Third  Basis  Function 

3 

<«> 

+  1 

N 

3 

-1 

N 

3 

(£,<pX<p,i) 

+1 

Y 

3 

(TJ.TJ) 

+  1 

N 

3 

(V>9)(9>V) 

-1 

Y 

3 

( <p,(p ) 

+1 

N 
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APPENDIX  I  -  COMPUTATION  OF  T-MATRIX  FOR  PROLATE  SPHEROID 


This  appendix  describes  the  method  for  performing  the  surface  integrals  necessary  to 
compute  the  spheroidal  T-matrix  for  a  prolate  spheroid.  Since  the  spheroid  is  both  symmetric 
under  the  transformation  Tj  -77  and  (p->(p+a,  the  T-matrix  is  diagonal  in  parity  <7  and 
azimuthal  order  m.  In  addition,  only  terms  for  which  the  difference  (/-/  *)  is  even  have  non-zero 
elements  for  the  T-matrix.  These  symmetries  allow  us  to  reduce  the  two-dimensional  surface 
integrals  to  a  one-dimensional  integral  over  the  variable  0  >  r\  >  1 .  Due  to  the  above  symmetries, 
the  T-matrix  is  of  the  following  form: 


Define  the  following  basis  functions  by  replacing  the  trigonometric  functions  in  the  basis 
functions  by  the  exponential  functions  exp [±im(p\ ,  respectively. 

exp[±/m^>]  =  -'(Sweml  ±  iS\ml) 

4L  exp [±inup]  =  -j==  (Aweml  ±  iAwoml ) 

v±m,  exp[±/m^]  =  -j==  (VZ  ±  i  Vol, ) 

4n,i  exp[±/m^>]  =  -j=  (AZ  ±  iAZ ) 

VL  expt+^3  =  ±iKm,) 

C,  exp [±im(p]  =  -^(t(vZ)±it(yZ)) 
exp[±/w^]  =  -j= L=  (t(VZ )  ±  it(VZ )) 

The  functions  Sw ,  Aw ,  Vw ,  and  tw ,  are  the  scalar  harmonic,  vector  harmonic,  vector  basis 
function,  and  traction  vector  in  the  external  fluid,  respectively.  The  functions  AT,  V1,  and  tT  are 
the  vector  harmonics,  vector  basis  functions,  and  traction  vectors  in  the  interior,  respectively. 
Depending  upon  whether  a  given  function  F±ml  contains  an  even  or  odd  number  of  derivatives 
with  respect  to  the  angle  (p  the  function  will  be  proportional  to  either  “1”  or  ±i . 
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Define  the  following  notation: 

p  [+1,  even  azimuthal  parity 
m/ 1  ±i,  odd  azimuthal  parity 


Define  the  following  functions  based  on  their  azimuthal  parity: 


zw  __  otr 
^ml  ~  ° ±ml 


A),=(4,,V± 

Ah=A,h 

Ah=(A,h 


(•£),=  (4*  V±' 

(^w)f  ~(Xnlh 


Ah=^Lh'±‘ 
Ah =(yL)„/±i 
A,h=A.,\ 
Ah  =(Oi 
Ah  =<0, 
Ah=A„,),'±i 
Ah=A„,h 
Ah  =(0, 
Ah=(Ah/:ti 
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:('L,v±/ 

:(tLi)n/±i 

'  (^±m/  )«7 

‘(O, 

1  (^±m/  );? 

1  (^±m/)» 7 

(£),  =(*L)f/:fcl 


<&)* 

<&), 

<S), 

(5)f : 
(S), 


The  computation  of  the  T-matrix  requires  the  evaluation  of  the  following  surface  integrals 
over  the  surface  of  the  target: 


e„v  =  -'-Jr  /<«{(»•  X»  •  4>» 


M  .  = 

n,n 


\dA{(n.V?)Sw,) 


P  . 


Jd4{(n»  ReVrn)SIFfr.} 


R,,,.,  =  famer,.)' 4J 

s 

In  the  above  integrals  one  adopts  the  convention  that  the  left  hand  term  in  the  integral  has  +m 
and  the  right  hand  term  has  -m  for  the  azimuthal  angle  dependence  of  the  integrand.  The  above 
integrals  are  diagonal  in  m,  and  they  can  be  replaced  by  the  following  one-dimensional  integrals: 

i) 

h,= -w-n1) 


=  h  h  n*  —  h~h  n * 

7/  V  £  g  <p  7] 


dr} 
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QT^-Uti-pr  \dr) 


-1 


M’,,  =  -2Ki  ]^A.(0{(51.) 


=  2*j<fljA,(ReO<(0 


=  2«r  }*;A„(ReO*(-0 

-I 

Using  the  polar  parity  and  the  azimuthal  parity  properties  of  the  integrands,  the  non-zero 
components  of  the  above  integrals  in  the  case  /-/  ’  is  even  are  given  by  the  following  integrals: 
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ev = ^4r]dn\az 
A  0 


=  4*J*,A„  (Re  VZ,)s(S:r) 
0 


P~w.=Az\dnA,^V^(Sl.) 

0 


^  MM'  “  4jr|./)7A„[+(Re?J,),A.)i+(Rei2,)«(i)»] 
0 


Rm2,,2r=4x  J [+(Re £ +  (Re t2m, % (A2mV ), ] 
0 

=4*  /^.(ReC^A), 

0 


Rm2/  3/.  =4^r  J  Jt7A„  (Re  Pm,  (A^ 
0 


=  4/}rf^a,[+(Re2,),(^.),  +(Re51)s(4L.),] 
0 
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The  non-zero  components  of  the  integrals  in  the  case  /-/’is  odd  are  given  by  the  following 
integrals: 


+4*i ]<*7A„(Re  V'„ )( (S'  ) 

0 


+1 

Rmu,ir  =  +A7ti  ldr>\[+(ReIlih,A)n 
0 


RmU,V=+iAj:  j dV^n  (Re  t'ml  )|  (^mr )| 

0 

In  the  case  (/-/  7  is  even,  the  components  of  the  R  and  P  matrices  are  real.  In  the  case  (/-/  0  is  odd 
the  R  and  P  matrices  are  imaginary.  The  R  and  P  matrices  can  be  transformed  into  real  matrices 
by  making  the  following  similarity  transformation  on  the  vector  indices  of  the  R,  P,  and  Q 
matrices: 


Sn.=-i5l'8)+8l'8]+8l'8] 


S-'n'=S'n.=+i8l'8)  +  8l'8]+8l'8] 


K,rr  =  I 


r",r" 


nm  _  Y  C  Pm 
rTU'  ~  Zj  °T\T'rT'U' 


In  the  case  (/-/’)  is  even,  the  components  of  the  Q,  R,  and  P  matrices  are  left  invariant  by  this 
transformation,  whereas  in  the  case  (/-/  ’)  is  odd  they  are  transformed  to  the  following  real 
matrices: 
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\i,r 


+1 


=  -(4*Jrf(j4,(ReKj,)<(S') 
0 


Rmuv.  =  +4k  ^dif&.n  (Re  Vm,  )^(A*r  ), 

0 

To  convert  the  above  integrals  into  matrices  one  adopts  the  following  convention  for  the 
index  of  the  pair  tI  of  vector  and  order  indices  in  the  case  m— 0,  where  negative  indices  are 
ignored: 

n(r  =  \,l)  =  3{l-m)-2 
n(r  =  2,1)  =  3(1 -m)— l 
n(T  =  3,l)  =  3(l-m) 

In  the  case  m*  0 ,  the  following  index  notation  for  the  pair  tl  of  vector  and  order  indices  is 
adopted: 

n(z  =  1,  /)  =  3(1  -  m)  +  2 
n(z  =  2,l)  =  3(l-m)  +  \ 
n(T  =  3,l)  =  3(l-m) 

Here  one  uses  the  C-style  index  notation,  where  the  first  element  in  a  vector  starts  at  the  index  0. 
The  above  index  convention  is  used  to  project  out  the  zero  components  of  the  vector  basis 
functions  in  the  case  l=m=0.  For  example,  the  square  R-matrix  has  dimension  (3(/max  -m)+ 1) 

in  the  case  m- 0,  and  (3(/max  -m)  +  3)  otherwise,  where  /max  is  the  maximum  value  for  the  order  /. 

The  T-matrix  is  defined  by  the  following  equation  in  terms  of  the  above  matrices: 

T  =  -(ReQR-'P-ReM)(QR-'P-My' 
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APPENDIX  J  -  COMPUTATION  OF  T-MATRIX  FOR  FINITE  CYLINDER 


This  appendix  describes  the  computation  of  the  spheroidal  T-matrix  for  a  finite  cylinder 
with  hemi-spherical  end  caps. 

First,  the  author  begins  with  a  discussion  of  the  geometry  of  the  surface  of  the  cylinder. 
Let  L  denote  the  length  of  the  cylinder  and  D  its  diameter,  where  the  aspect  ratio  %  =  (LID)  is 
greater  than  unity.  The  semi-focal  distance  of  the  spheroidal  coordinate  system  is  defined  by  the 
following  equation  in  terms  of  the  length  and  aspect  ratio  of  the  cylinder: 

/=(i/2)^E> 

X 

Define  the  quantity  Lcy,  -L-D  as  the  length  of  the  cylindrical  portion  of  the  target 
excluding  the  end  caps. 

In  the  case  —Lt  /  2  <  z  <  +Lcyl  /  2 ,  the  normal  to  the  surface  is  given  by  the  following  equation: 


(ncyi)f  -+£ 


l-?72 

f-V2 


( ncyl)r)=-V . 


cf-1 

f-t?2 


(«c A=° 

The  spheroidal  coordinates  of  the  point  on  the  surface  in  the  <p  =  0  plane  are  given  by  the 
following  equations: 
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p  =  plf  =  Dlfl  2 
z  =  z!  f 


■n=z/% 

To  describe  the  surface  of  the  target  at  its  end  caps  one  uses  the  polar  angle 
-n  /  2  <  <  +Jt  /  2  to  parameterize  the  surface  at  its  end  caps.  The  coordinates  of  a  point  on  the 

end  caps  in  the  <p  =  0  plane  are  given  by  the  following  relationships: 

p  =  (D/2/)sin(tf) 
z  =  z/f  =  ±((D  /  If)  cos(tf)  +  Lcyl  /  2) 


ri  =  zl% 

The  normal  to  the  surface  at  its  end  caps  in  the  (p  =  0  plane  is  given  by  the  following  equations: 
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Z  =  ffr 

nx  =  sin($) 
nz  ~  cos(i}) 


n{  =  nxdsx+n,dfz  =  sin(tf)  +  77  cos(tf)] 

nn  =  nxdnx + n2dnz  =  f[-jj  sin(*?)  +  €  cos(0)] 

np=0 


tl;  =  [+^\-7J2  sin(7?)  +  7]'I¥~- 1  cos(i})] / yj<^2 -T]2 
Hf,  =  [~W^2_1  sin(^)  +  COS(tf)] / V£2-772 


The  component  7^  of  the  normal  has  even  parity  under  the  transformation  77  -4  -77 , 
whereas  the  component  has  odd  parity.  In  the  case  of  the  normal  for  the  end  caps,  the  angle 
z?  becomes  t?  — >  d+n  under  the  transformation  77  -4 -77 . 


Since  the  finite  cylinder  is  both  symmetric  under  the  transformation  77-4-77  and 
(p  — >  (p+a ,  the  T-matrix  is  diagonal  in  parity  <7  and  azimuthal  order  m.  In  addition,  only  terms 
for  which  the  difference  (/-/’)  is  even  have  non-zero  elements  for  the  T-matrix.  These 
symmetries  allow  us  to  reduce  the  two-dimensional  surface  integrals  to  a  one-dimensional 
integral  over  the  variable  0  >  77  >  1 .  Due  to  the  above  symmetries,  the  T-matrix  is  of  the 

following  form: 


np  _  Off'  Cm'rrim 

Ioml,o'mT  —  °a  °m  1l,l’ 


Define  the  following  basis  functions  by  replacing  the  trigonometric  functions  in  the  basis 
functions  by  the  exponential  functions  exp[±w»^>] ,  respectively: 
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Sw±m,  exp [±im<p]  =  -j=^(SH' eml  ±iSw (mj) 

Kmi  exp [±im(p\  =  ± iA*ml) 

VL  exp[±imp]  =  -j===(V"ml  ±  iVZ,) 

Ar±ml  exp [±im<p]  =  -j==(Alml  ± iAToml) 

VL  exp [±im(p]  =  -j=(VeTml  ± iVzoml) 
tlm,  exp[±im<p]  =  -j==  (t(Veh )  ±  it(VZi )) 
tT±m i  exp [±im(p]  =  WL )  ±  *'(C/ )) 

The  functions  S**',  Aw ,  Vw ,  and  tw ,  are  the  scalar  harmonic,  vector  harmonic,  vector 
basis  function,  and  traction  vector  in  the  external  fluid,  respectively.  The  functions 
Ar ,  VT ,  and  f  are  the  vector  harmonics,  vector  basis  functions,  and  traction  vectors  in  the 
interior,  respectively.  Depending  upon  whether  a  given  function  F±ml  contains  an  even  or  odd 
number  of  derivatives  with  respect  to  the  angle  (p ,  the  function  will  be  proportional  to  either  “1” 
or  ±i. 


Define  the  following  notation: 

~  f+1,  even  azimuthal  parity 

±m/  m/ j±/5  odd  azimuthal  parity 

Define  the  following  functions  based  on  their  azimuthal  parity: 


rW  _  oW' 
°m/  —  ° ±ml 


(4 L,),=(4.,V±<’ 
A\=(4„), 

(4»,){=(-Os 
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(v:,\={viml\i±i 

(vL\=(vL,X 

<yi)t=(rL)t 

(riX=<rLx 

<yih=<yL)t 

(vl,\=<yLx 


(*ml\  ~  (*±ml)e  ~  ! 

&), =(>D,'±- 

(fmlX  ~  (t±m/)p 
=(4m/)| 
(XlX  =  (*±mlX 
(%),=&),  /±l 
(*mt  = 

(XlX  ~  (t±ml \ 
(<»,)»  =«LV±i 


The  computation  of  the  T-matrix  requires  the  evaluation  of  the  following  surface  integrals 
over  the  surface  of  the  target: 
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Q,.r,  =-i-^w  \dA’(n*i(V“')(H‘Ar._.)} 

K,=-i-pr  /<««»' <)*’.■) 

A  s 

P„.,=  jdA«n'ReV,,)S\.} 

s 

R„„=\dAmeV,J*A,.„.) 

s 

In  the  above  integrals  one  adopts  the  convention  that  the  left  hand  term  in  the  integral  has  +m 
and  the  right  hand  term  has  —m  for  the  azimuthal  angle  dependence  of  the  integrand.  The  above 
integrals  are  diagonal  in  m,  and  they  can  be  replaced  by  the  following  one-dimensional  integrals: 


Pm„,  =  2*  \dn A.KReOrt  +<R*0,»i,K0 


=2*  J<*7A„(R  et;j.(4) 

“I 

Using  the  polar  parity  and  the  azimuthal  parity  properties  of  the  integrands,  the  non-zero 
components  of  the  above  integrals  in  the  case  /-/  ’  is  even  are  given  by  the  following  integrals: 
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P\,r  =  471  jdn\  [+(Re  K2, )(  »(  +  (Re  Vi )  „  it,  X& ) 
0 


P’w  =  4ir  Jdl7A,[+(Re^B|  +(Refi)4»,)(^.) 
0 


0 


Rm2,,2r  =  ^  \dTlK [+(Re i2m, )„  A  )*  +  (Re Pm, )^A2m,,)^] 
0 


*"w=^]rf>;A.(ReC){(i)| 

0 


+1 

Rm 2i,2i'  =  \dV\  (Re  4/)|(4mr  )| 

0 


=  4*  Jrf7A,[+(Re?i),(i),  +(Re»i), 
0 
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The  non-zero  components  of  the  integrals  in  the  case  /-/  ’  is  odd  are  given  by  the  following 
integrals: 


0"/f„.  =  -4 ^n'(-o4r  prtjA.WE +  (C 
^  0 


Pmur  =+AKi\dvK[^Vl,)^  +(ReFi)^](S.) 
0 


=+4»;  j^Ai+(Re?;,),(i)4  -(Re<w)#(i)»] 
0 


*"„.>/•  =  ^  }rf/?A,(Re?i)^  (i.)j 
0 

In  the  case  (/-/  0  is  even,  the  components  of  the  R  and  P  matrices  are  real.  In  the  case  (/-/  “)  is  odd 
the  R  and  P  matrices  are  imaginary.  The  R  and  P  matrices  can  be  transformed  into  real  matrices 
by  making  the  following  similarity  transformation  on  the  vector  indices  of  the  R,  P ,  and  Q 
matrices: 


S7,=-iFr'%+8;'%+8l'% 


S-\,  =  5+rr,  =  +iSzT'%  +  8TZ'%  +  8Tr% 


—  V  c 

KtIjT  “  Zj  ^r',r"/V7,r’V°  rV 


r",r" 


j 


P;,='LSr.,-P?,.r 

X* 


x' 
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In  the  case  (/-/ ")  is  even,  the  components  of  the  Q,  R,  and  P  matrices  are  left  invariant  by  this 
transformation,  whereas  in  the  case  (/-/’)  is  odd  they  are  transformed  to  the  following  matrices: 

fi-w=-4*i'4r 

^  0 

=-(4*  +(Re^)4«,)(S;.) 

0 

+1 

R  U2r  =  +4;r 

o 


0 

To  convert  the  above  integrals  into  matrices  one  adopts  the  following  convention  for  the 
index  of  the  pair  tI  of  vector  and  order  indices  in  the  case  m- 0,  where  negative  indices  are 
ignored: 

n(T  =  1,  /)  =  3(/  -m)-2 
n(r  =  2,/)  =  3(/-m)-l 
n(r  =  3,l)  =  3(l-m) 

In  the  case  m  *  0 ,  one  adopts  the  following  index  notation  for  the  pair  tl  of  vector  and  order 
indices. 

n{t  =  1,  /)  =  3(7  -m)  +  2 
n(t  =  2,  /)  =  3(1  -m)  +  \ 
n(r  =  3,l)  =  3(l  -m) 

Here  one  uses  the  C-style  index  notation,  where  the  first  element  in  a  vector  starts  at  the  index  0. 
The  above  index  convention  is  used  to  project  out  the  zero  components  of  the  vector  basis 
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functions  in  the  case  l=m~0.  For  example,  the  square  R-matrix  has  dimension  (3(/max  -m)  + 1) 
in  the  case  0,  and  (3(/max  -m)  +  3)  otherwise,  where  /max  is  the  maximum  value  for  the  order  /. 

The  T-matrix  is  defined  by  the  following  equation  in  terms  of  the  above  matrices. 

T  =  -(Re  QR-'P  -  Re  M)(QR~lP-  M)'1 
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